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ABSTRACT 

Aims. We present the size, shape and flux densities at millimeter continuum wavelengths, based on ALMA science veriflcation 
observations in Band 3 (~ 94.6 GHz) and Band 6 (~ 228.7 GHz), from the binary Mira A (o Ceti) and Mira B. 

Methods. The Mira AB system has been observed with ALMA at a spatial resolution of down to ~ 25 mas. The extended atmosphere 
of Mira A and the wind around Mira B sources are resolved and we derive the size of Mira A and of the ionized region around Mira 
B. The spectral indices within Band 3 (between 89-100 GHz) and between Band 3 and Band 6 are also derived. 

Results. The spectral index of Mira A is found to change from 1.71+0.05 within Band 3 to 1.54 + 0.04 between Band 3 and 6. The 
spectral index of Mira B is 1.3+0.2 in Band 3, in good agreement with measurements at longer wavelengths. However it rises to 
1.72 + 0.11 between the bands. For the first time the extended atmosphere of a star is resolved at these frequencies and for Mira A 
the diameter is ~ 3.8 x 3.2 AU in Band 3 (with brightness temperature ~ 5300 K) and ~ 4.0 x 3.6 AU in Band 6 (T^ ~ 2500 K). 
Additionally, a bright hotspot ~ 0.4 AU and with ~ 10000 K is found on the stellar disc of Mira A. The size of the ionized region 
around the accretion disk of Mira B is found to be ~ 2.4 AU. 

Conclusions. The emission around Mira B is consistent with that from a partially ionized wind of gravitationally bound material from 
Mira A close to the accretion disk of Mira B. The Mira A atmosphere does not fully match predictions, with brightness temperatures 
in Band 3 significantly higher than expected, potentially due to shock heating. The hotspot is likely due to magnetic activity and could 
be related to the previously observed X-ray flare of Mira A. 
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1. Introduction 


Mira i s the closest symbiotic binary star (at 92 pc. Ivan LeeuwenI 
(l2007l) ) consisting of a regularly pulsating, mass-losing Asymp¬ 
totic Giant Branch (AGB) primary, Mira A (o Get), and a com¬ 
panion, Mira B (VZ Get), believed to be a white dwarf. Resolved 
UV observations of th e two components were first published in 
Karovska et al.l (Il997l) and revealed material flowing from Mira 
A to Mira B. 

A soft X-ray outb urst was detected from Mira A in 2003 
jKarovska et al.ll200^ . Since then, the system has been closely 
monitored at different wavelengths to further understand the 
cause and effects of the X-ray outburst. The OH maser emis¬ 
sion is also varying a nd the most recen tly observed OH maser 
flare occured in 2009 (lEtoka et al.ll2010l) . Although the OH flar¬ 
ing can be correlated with variations in the optical lightcurve and 
the UV line emission occuring on the same time scales, a con¬ 
sistent model explaining the transient phenomena detected by 
different probes is still missing. 

The submm emission from the Mira AB system around the 
CO(J =3^2) line at 345 GH z was observed with ALMA in Cy¬ 
cle 1 (iRamstedt et al.|[20T^ . The CO maps show that the cir- 
cumstellar gas has been shaped by several different dynamical 
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processes during the evolution of the star. In the continuum emis¬ 
sion centered on 338 GHz, the binary pair is marginally resolved. 
In this paper, we present new ALMA continuum observations of 
the system taken as part of the ALMA long-baseline science ver¬ 
ification campaign (ALMA partnership, 2015). 


2. Observations and analysis 

Mira was observed on 17 and 25 October 2014 (Band 3) 
and 29 October and 1 November 2014 (Band 6). Here we 
present results from the dual polarization continuum spectral 
windows. In Band 3, three continuum spectral windows were 
observed with 128 channels over 2 GHz bandwidth centered 
on 89.176 GHz, 91.234 GHz and 99.192 GHz. In Band 6, one 
continuum 2 GHz spectral window with 128 channels was cen¬ 
tered at 228.558 GHz. We used the data products provided by 
the ALMA observatory except that the two epochs for each ob¬ 
serving band were treated separately as the sources are strongly 
detected in each epoch, allowing for the study of structural and 
flux changes. Self-calibration was performed on the continuum 
spectral windows that were averaged from 128 to 8 channels 
across the 2 GHz bandwidth. Flux density calibration was per¬ 
formed using a number of different quasars for the two bands and 
epochs. The calibrators are listed in Appendix I aI Based on the 


Article number, page 1 of|5] 


















A&A proofs: manuscript no. Vlemmings_Mira 



0.6 0.4 0.2 0 - 0.2 100 50 0 -50 100 50 0 -50 

East offset [arcsec] East offset [mas] East offset [mas] 


Fig. 1. Continuum image of the Mira AB 
binary at Band 6 (left). Residuals in Band 
6 at the location of Mira A after subtract¬ 
ing the best fit disc model (middle) and the 
disc model including a compact gaussian 
hotspot (right) from the visibilities. The 
dashed ellipse indicates the size of the fit¬ 
ted stellar disc. Contours in all three plots 
are drawn at 3,30,90,270, 810cr with cr = 
140 pJy beam"^ The beam size is shown in 
black in the bottom left corner of all three 
figures. 


uncertainty of the catalog values due to quasar variability and the 
flux densities measured for the phase and bandpass calibrators 
in the different data sets, we conservatively estimate an absolute 
flux density calibration uncertainty of ~ 5% in both observing 
bands. 

Most of the data analysis was performed directly on the uv- 
data in order to avoid image deconvolution artifacts, but im¬ 
ages were also created using the CASA task clean with Briggs 
0.5 weighing of the visibilities. This resulted in a beam of 
72 X 64 mas and a position angle of 78.8° on October 17 and 
68 X 57 mas at a position angle of 73.3° on October 25 for 
Band 3. For Band 6 both days reached the same resolution of 
34 X 24 mas with a position angle of 22.4° on October 29 and 
16.0° on November 1. The rms sensitivity was ~ 40 yuJy beam“^ 
in both epochs of the Band 3 and ~ 130 yuJy beam“^ in both 
epochs of the Band 6 observations. The images in Band 3, with 
> 10% fractional bandwidth coverage, were made using a the 
multi-frequency synthesis option in clean with two Taylor coef¬ 
ficients to model the frequency dependence. This allowed us to 
determine the spectral index at the ~ 94.2 GHz representative 
frequency. 

T o analyze the t/v-data dir ectly, we used the UVMULTIFIT 
code (iMarti-Vidal et al .1120141) . We And that the observations of 
October 25 in Band 3 are likely affected by a marginally re¬ 
solved quasar flux calibrator, producing systematically smaller 
sizes, and do not consider this epoch in our analysis. 


3. Results 

We present a representative image at the two frequencies in 
Fig- [B (left, band 6) and Fig. [2] (left, band 3). The relative as¬ 
trometry of the observations is discussed in Appendix |B] In or¬ 
der to determine the size and shape of both Mira A and Mira B, 
we performed i/v-fltting. For both sources we attempted fits to 
unresolved delta functions, elliptical and circular Gaussians and 
elliptical and circular discs. The best fit parameters for Mira A 
and B at each epoch and observing band are presented in Ta¬ 
ble [T] As indicated in Fig. [2| (right), the emission from Mira B 
is clearly extended and the best fits were obtained using a cir¬ 
cular Gaussian with a spectral index of a = 1.3. Within the fit¬ 
ted errors and estimated absolute flux density uncertainties, the 
flux density of Mira B is roughly constant at the two epochs in 
both bands, although in Band 3 there might be a hint of vari¬ 
ability at the 5 - 10% level. Between Band 3 and Band 6, the 
spectral index is 1.72 + 0.11, taking into account all flux den¬ 
sity uncertainties. Ignoring the observations of October 25, the 
size is consistent between the different frequencies with a full- 
width half-maximum (fwhm) of 26 mas, which corresponds 
to ~ 2.4 AU. 

For Mira A, the best fits to the Band 3 data were produced 
with an elliptical disc with a ^ 1.7. The flux density is nearly 


constant during both days at each observing frequency. The spec¬ 
tral index between Band 3 and Band 6 however, is shallower 
with a = 1.54 ± 0.04. The disc is clearly elongated, with a major 
axis of 41.8 mas (3.8 AU), axis ratio of ~ 1.2 and position angle 
~ 54°. However, in Band 6, we found as indicated in Fig.[T] a sig¬ 
nificant component remaining after subtracting the best fit disc 
model. This component, consistently at both epochs, could be fit 
by an additional compact gaussian offset by ~ 3 mas from the 
disc center. Fitting both disc and gaussian simultaneously signi- 
cantly reduced the residuals in the image. The size of the disc 
is larger at Band 6 compared to Band 3, with an average size of 
43.32 mas (3.99 AU) and an axis ratio of ~ 1.1. 


4. Discussion 


4.1. Mira A 

4.1.1. Size and shape 

Mira A has long been known to have an asymmetric stelly disc 


(e.g.lKarovska et al.ll99lllHaniff et al.l 19921: 1 Wilson et al.ll992t 


lOuirrenbach et al.l[l992l) . The position angle and axis ratio how¬ 

ever, has been shown to be quite variable across the stellar phase. 
Most observations taken around maximum phase are found to 
have a position angle of the major axis between ~ 105 - 160°. 
Closer to the phas e of the ALMA observ ations (phase ~ 0.36) 
the observations of iKarovska et al.l (Il99ll) (at phase 0.43) And a 
position angle of ~ 30°. The elongation itself varies with wave¬ 
length and phase from none to ~ 20%, consistent with our obser¬ 
vations. The wavelength dependence of the elongation is likely 
due to opacity effects. The apparent variation of position angle 
with stellar phase could indicate that stable non-radial pulsations 
shape the extended atmosphere. 

The diameter of Mira A also varies with phase by about 
10% and was measur ed to be ~ 32 mas at 2 pm wavelength 
(I Woodruff et al.ll2004b . This corresponds to a Rosseland radius 
of ~ 1.5 AU. The radius derived from the major axis we mea¬ 
sure at Band 3 would thus corresponds to ~ 1.3 while at 
Band 6 it corresponds to ~ 1.35 R^. Previous observations at 
43 GHz resolved the radio photosphere of Mira A and derived a 
uniform disc diameter of ~ 52 mas, corresponding to ~ 1.6 R^ 
iReid & Mentenl(l2007l) . The fact that the Band 6 observations ap¬ 
pear to probe a larger radius than the Band 3 observations likely 
requires different opacity sources throughout the atmosphere. 


4.1.2. Emission process 

iReid & MentenI (1 19971) analyze different centimeter-wavelength 
continuum emission processes and conclude that the detected 
emission at 8.5-22 GHz originates in a radio photosphere reach¬ 
ing unity optical depth at around ~2R^, predicting spectral in- 
dex of = 1.86 at centimeter wavelengths. iMatthews & Karovskal 
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Fig. 2. Continuum image of the Mira AB 
binary at Band 3 (left). Residuals are shown 
for Band 6 at the location of Mira B af¬ 
ter subtracting the best fit delta function 
(middle) and extended Gaussian (right). 
The dashed circle indicates the fwhm of 
the fitted Gaussian. Contours in all three 
plots are drawn at 3,30,90,270,810cr with 
cr = 40 yuJy beam“^ in Band 3 and cr = 
140 beam“^ in Band 6. The beam size 
is shown in the bottom left of all three fig¬ 
ures. 


Table 1. wv-fitting results 


epoch 

y 

[GHz] 

shape 

[mJy] 

major axis / fwhm 
[mas] 

axis ratio 
major/minor 

position angle 

n 

spectral index 

Mira A 

17 Oct 2014 

94.2 

Disc 

35.03+0.04 

41.8+0.4 

1.20+0.01 

54+2 

1.73+0.09 

25 Oct 2014^ 

94.2 

Disc 

34.52+0.04 

38.6+0.4 

1.28+0.01 

54+1 

1.70+0.04 

29 Oct 2014 

228.67 

Disc 

137.8+0.2 

43.28+0.07 

1.13+0.02 

51.0+0.5 



228.67 

Gaussian 

10.13+0.07 

4.6+0.5 

1.0 



01 Nov 2014 

228.67 

Disc 

140.0+0.2 

43.36+0.06 

1.12+0.02 

50.8+0.6 



228.67 

Gaussian 

8.98+0.07 

4.7+0.5 

1.0 



Mira B 

17 Oct 2014 

89.06 

Gaussian 

2.50+0.04 

24+2 

1.0 


1.3+0.2 

25 Oct 2014^ 

89.06 

Gaussian 

2.25+0.04 

18+2 

1.0 


1.3+0.3 

29 Oct 2014 

228.67 

Gaussian 

12.16+0.09 

25.7+0.3 

1.0 



01 Nov 2014 

228.67 

Gaussian 

11.98+0.09 

26.2+0.3 

1.0 




^ The size of both Mira A and B are likely underestimated at the second epoch due to a marginally extended flux calibrator J2258-2758 


(l2006h published resolved images of the Mira system at 8.5- 
43.3 GHz. The system was monitored with the VLA over ^^80% 
of the pulsational cycle and showed flux density variability 
<30% below 22 GHz. T he lower frequency da ta is consis¬ 
tent with the model from iReid & Mente^ (1 19971) . At 43 GHz 
the flux density ap pears to show larger v ariability based on 
the observations by iReid & MentenI (l2007l) who found a flux 
density almost a factor of t wo higher than that reported in 
[Matthews & Karov^ (l2006l) . The VLA values, the ALMA 
values from this wo rk and the flux density at 345 GHz 
(iRamstedt et al.ll20T4l) are presented in Fig.O The most striking 
result is that the spectral index between the ALMA Band 3 and 
Band 6 observations (a = 1.54+0.04) appears signiflcantly more 
shallow than that within Band 3 itself (a = 1.70 ± 0.05). The 
shallower slope agrees with the observations in ALMA Band 7 
though the steeper slope would be required to fit the radio ob¬ 
servations. Both spectral indexes are not quite in agreement with 
the model prediction bv iReid & MentenI (1 19971) . indicating that 
at millimeter wavelengths the sources of opacity change. Our ob¬ 
servations yield signiflcantly higher brightness temperature than 
expected from the VLA observations. For a uniform disc we And 
Tb ~ 5300 K at Band 3 (~ 1.3 R.) and Tj, ~ 2500 K at Band 6 
1.35 R,), compared to then ~ 1650 K at 43 GHz(~ 1.6 R,). 
Only part of this difference could be explained by variability and 
especially the higher temperatures in Band 3 might require shock 
heating of the atmosphere close to the stellar surface. 

4.1.3. Stellar Activity 

As shown in Figdl the Band 6 data require the presence of a 
strong compact component with a flux density of ~ 10 mJy off¬ 
set by ~ 3 mas from the stellar disc center. A gaussian compo¬ 


nent of ~ 4.7 mas (~ 0.4 AU) produces a good fit, though cannot 
rule out a cluster of delta functions. We have thus detected the 
presence of a significant hotspot, or a compact cluster of spots, 
on the Mira A stellar disc at ~ 1.3 mm wavelength. We have 
investigated if the same spot could be detected in the Band 3 ob¬ 
servations at ~ 3 mm. Fits including a gaussian component up to 
~ 2.5 mJy produced equally good results, but the worse angular 
resolution did not allow us to confidently distinguish between 
a fit of a stellar disc plus compact component or a stellar disc 
alone. 

The flux density measurements allow us to determine the 
brightness temperature. With a size of ~ 4.7 mas, we And 
Tb ~ 10000 K, above the brightness temperature of ~ 2500 K 
measured for the stellar disc in Band 6. The upper limit in Band 
3, assuming a 4.7 mas area, gives Tb < 17500 K. Bright hotspots 
at millimeter wavelenghts could be caused by shock heating due 
to pulsations or convection. However, the high brightness tem¬ 
perature of the hotspots at ~ 1.3 R^, is more readily explained 
by magnetic activity as seen on our Sun. Similar magnetic flares 
were also suggested to be the cause of t he soft X-ray outburs t 
observed on Mira A in December 2003 (iKarovska et al.|[200^ 
but this is the first direct detection of such magnetic activity 
of an AGB star in the long wavelength regime. The observed 
hotspot could be related with the strongly polarized elongated 
SiO n iaser ejections that a ppear aligned w i th a ra dial magnetic 
field (ICotton et al.l l2006l) . IKarovska et al.l (l2005l) crudely esti¬ 
mated the lifetime of such magnetic flares to be of order several 
weeks. We notice a significant difference of almost 15% within 
2 days. This would be consistent with the lifetime estimate of 
a few weeks, but could also indicate significant short term vari¬ 
ability. 
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4.2. Mira B 

The accreting companion Mira B has a time-varying accretion- 
driven wind (M ~ 5 X 10~^^ - - 10" ^ M© yr“^, Voo ~ 250 - 
450kms“\ IWood & Karovskal l2006h and recently ALMA 
showed the impact of this fast and tenuous wind on the much 
slower and denser outflow (M ~ 10“^ M© yr“\ Vc» ~ 10 km s“^) 
of Mira A (iRamstedt et al.ll20T4l) . Even if fully ionized, the ex¬ 
pected flux density of Mira B ’s wind would be at least 4 orders of 
magnitude smaller than that we observe at 89 and 229 GHz, and 
would also be optically thin (i.e., Sy oc We therefore rule 

out any significant contribution from Mira B’s wind at ALMA 
frequencies. The steep spectral index of Mira B (tr ~ 1.3 - 1.7) 
also rules out any significant non-thermal component, while the 
symmetric nature of the resolved emission makes it unlikely that 
the emission is being produced from a collision of the two winds. 
Our derived diameter for Mira B is ~ 2.4 AU which is 2 orders 
of magnitude greater than the si ze of the accretion disk detected 
bv iReimers & Cassatellal (Il985l) and so we do not directly detect 
the accretion disk at ALMA frequencies. 

[Matthews & Karovskal (l2006l) treated the centimeter emis¬ 
sion from Mira B to emanate from an ionized cavity in Mira 
As wind, formed by the UV radiation field of the hot accre¬ 
tion disk. They approximated a diameter of 1 - 5 AU for Mira 
B, which is consi s tent w ith our findings. However, following 
iTavlor & SeaquistI (Il984l) and assuming a gas temperatures of 
T = 10"^ K, the spectral energy distribution for such an ionized 
cavity should have turned over to a spectral index of a = -0.1 
well below our ALMA frequencies. Instead we find that a = 
1.3 + 0. 2 at 89 GHz, which is consisten t with the a = 1.18 + 0.28 
value of iMatthews & Karovskal (l2006l) at lower frequencies. The 
absence of a turnover in the spectral index up to 89 GHz implies 
that the gas density at the location of Mira B is at the very least 
an order of magnitude greater than the gas density of Mira As 
wind at that location. It is therefore likely that we are detect¬ 
ing a gravitationally bound partially ionized gas at centimeter 
and (sub-)millimeter wavelengths and not an ionized cavity of 
Mira As wind. We also find that the spectral index steepens to 
a = 1.72 + 0.11 between 89 and 229 GHz, which could be due 
to additional dust emission at a few tens of K. 


5. Conclusions 

With the longest ALMA baselines we have, for the first time, 
been able to resolve both the extended atmosphere of Mira A and 
the partially ionized region around Mira B at millimeter wave¬ 
lengths. The partially ionized region of ~ 2.4 AU around Mira 
B requires a density at least an order of magnitude higher than 
expected from the Mira A wind. This supports the suggestion 
that the region is made up from gravitationally captured material 
from the AGB wind that will be accreted on Mira B through its 
accretion disk. The elongation and wavelength dependent size 
of Mira A point to strong changes of opacity source through¬ 
out the extended envelope and possible non-radial pulsations. 
The brightness temperature at 3 mm is significa ntly higher than 
predic ted in the radio photosphere model from iReid & MentenI 
(l1997h . At 1.3 mm wavelength we also find a hotspot with an 
area of ~ 8% of the stellar disc. The hotspot has ~ 10000 K. 
This indicates an origin likely connected to magnetic activity and 
potentially related to the process responsible for the previously 
observed X-ray outburst and other recorded transient phenom¬ 
ena. The stellar activity could also explain the hot layers of the 
atmosphere probed by the 3 mm observations. 



V (GHz) 


Fig. 3. The Mira A and B spectral energy distribution from the radio to 
the sub-millimeter frequency range as probed by ALMA and the VLA. 
The round (red) symbols denote Mira A and the squares (blue) Mira 
B. The grey symbols are literature referred in the text. The errors on 
the VLA data points only contain the formal flux density uncertainties, 
while the ALMA error bars also contain the absolute uncertainties (5%, 
and 10% for Band 3/6, and 7 respectively). The dotted lines indicate the 
iReid & MentenI jl997l) model for Mira A and the fit to the VLA data 
for Mira B. The dashed line at the Band 3 points is the spectral index 
derived within the Band and the green dashed line is the fit to the ALMA 
Band 3 and 6 data. 
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Table A.l. Calibrator fluxes 


epoch 

calibrator 

reference freq. 
[GHz] 

flux density 

[jy] 

spectral index 

17 Oct 2014 

J0334-4008 

86.23 

1.66 

-0.71 

25 Oct 2014 

J2258-2758 

86.23 

1.21 

-0.73 

29 Oct 2014 

J0334-4008 

229.55 

0.83 

-0.70 

01 Nov 2014 

J0238+1636 

229.55 

1.39 

-0.20 


Appendix A: Flux density calibrators 


Several quasars were used for flux densi ty calibration at the dif¬ 
ferent epochs and bands. In Table jA.ll we present the calibra¬ 
tors, and their adopted flux densities and spectral indices from 
the ALMA calibrator catalog. 


Appendix B: Relative astrometry 

The binary pair Mira AB is fully resolved in the observations, 
with a representative image at the two frequencies presented 
in Fig. [T] (left) and Fig. [2] (left). Using wv-fltting we deter¬ 
mine the separation and position angle with milliarcsecond ac¬ 
curacy as indicated in Table ibh The average separation of 
0.472"corresponds to 43.4 AU at a distance of 92 pc. These val - 
ues flt well with the predictions made by iPrieur et ^ (l2002h . 
although the authors note that their derived binary orbit is still 
poorly constrained. It is interesting to note that the data already 
show a (not yet signiflcant) sign of a decrease in separation R 
over two weeks of observing at a rate of - 121 + 71 /ias day“^. 
Futher epochs of observations with ALMA will be thus be able 
to constrain the binary orbit to high precision. 

Table B.l. Relative astrometry 


epoch 

separation 

["] 

position angle 
[°] 

17 Oct 2014 

0.4722+0.0005 

98.79+0.07 

25 Oct 2014 

0.4721+0.0005 

98.59+0.06 

29 Oct 2014 

0.4719+0.0002 

98.64+0.02 

01 Nov 2014 

0.4709+0.0001 

98.53+0.02 
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